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APPROXIMATION 

\i  f\ 

DANS 

L’INTERPOLATION 


Soit  f(x)  une  fonction  qui  varie  constamment  dans  le  même 
sens  quand  x croît  de  i à N;  supposons-la  croissante  et  posons 
f(n)  = v„. 

Soient  V, , V2,  V3,  ...  des  multiples  de  io~*  ( k étant  un  entier 
positif)  respectivement  approchés  de  e,,  e2,  e3,  ...  à ^ to“*  près. 
Posons  pour  abréger, 

U — IO-'q  O = l-  M,  dn—  Vn+X  — Vni  D„  ==  V n-hi  — V*. 

Par  hypothèse, 

0)  | V;i  — en|  = 8 ; 

par  conséquent, 

(a)  | D„—  dn\ 4 i S — u. 

Donc  dn  — u — a , car  dn^>  o;  mais  Dn  est  un  multiple 
de  u , donc  D^o  ou  V„+1  > Vrt  ; la  suite  V,,  Y2,  V3-,  ...  est  /zo/i 
décroissante  : si  V„>>  VW',  il  s’ensuit  que  n S > 

Considérons  les  fonctions  fn{x)  et  Fre(#)  définies  par 

(3)  fn(x)—  p„H-  (x  — n)  dn=  p„+  (x  — n)(  v»), 

F»0)  = V«-t-  (#  -r-  n)  Dn—  — n)(VB+i- V»). 

Ces  fonctions  sont  croissantes  (sauf  que  F,j(&)  est  constante 

quand  D„  = o).  De  plus, 

f n (n)  — ( n)  = Vw,  f n{  n H-  x)  = vn-\-i , F^(  n t ) = \ n+\  ; 

( 4 ) fn  ( x ) fn  ( X ) = ( X X ) 0^, 

(5)  F„(#)—  F^f#')  = (-a?  — x')  Dn. 
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En  désignant  /«(a?)  — F«  (x)  par  gn(x),  on  a 

(6)  £nO)  = /«0)  — F„0) 

= (n-h  i — x)(vn  — \n)  -4-  (x  — n)(vn+x—  V»+1): 

donc,  d’après  (i), 


| gu{x)  | ^ 8 (|  n -+- 1 — |+  [ # — n I). 


Étudions  la  variation  de  | n + i — a:  | -j-  | a?  — n | quand  x croît 
de  n — i à /i  4-  2, 


1 

3 

X 

n — 1 

n n /i  + i 

2 

Tl  — H — 

2 

n H-  2 

| n + i — x 1 -h  \ x — n | 

3 

2 1 1 

2 

3 

donc 


(7) 

(8) 

(9) 


si 

si 


n % x £ n -h  i , 

1 3 

n <x<n-\ > 

2 ” 2 

n — i ^x^n  -h  2, 


I 8 n(æ)  | =.2  S, 

I £>'0)1^3  8. 


Si  l’on  admet  que,  pour  n — i < x «<  n -f-  2,  f(x)  — /,*(#)  soit 
négligeable,  on  pourra  écrire 

/O)  — P„0)  = fn(x)  — F„0)  .=  gn{x) 

et  alors  les  formules  (/),  (8),  (9)  donnent  une  limite  supérieure 
de  l’erreur  commise  en  remplaçant  f(x)  par  F„(a:). 

Si  n<fx<^n-\-  1 et  si  w est  une  valeur  approchée  de  Fn(a?) 
à û près,  on  a 

(10)  I/O)  — w | < 28  = u. 


Problème  inverse.  — b et  $ étant  des  nombres  donnés  (8  ]>  o), 
trouver  une  limite  inférieure  et  une  limite  supérieure  des 
nombres  x qui  vérifient  la  relation 

b — P C/O)  <b  h-  p ou  I/O)  — b | < [B. 


On  commence  par  chercher  des  entiers  ai,  n\  n tels  que 

b <V„+1, 

b — Ç>-6<  V„.+1,  V*.-! < 6 ■+■  p + 8 £ V*-. 
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On  en  conclut  d’abord  que  n'  <i  n 4-  i et  a <<  n \ ensuite 

o*4  V»'  + 8 < A - (3,  b -f-  g S Y*. ■ - S i *v. 

Donc  n'  <^x  <i  n"  et  l’on  s’en  tient  là  si  n'  < n — i ou  n"  n -f-  2 . 
Si  n — i < n'  et  n"£  n - f-  2,  on  a 

n — i < a?  < /i  h-  2 ; 

on  peut  alors  essayer  de  trouver  une  limite  inférieure  et  une  limite 
supérieure  de  x plus  resserrées  que  n'  et  n" , en  supposant  tou- 
jours f(x)  — fn  (x)  négligeable  pour  n — i <C  x <C  n,  -P  2.  Pour 
cela,  on  calcule  le  nombre  a défini  par 

F«(a)  = b. 

Gomme  Y nS  b <[  Vw+, , on  a 

n^a  <C  n - h i . 

fn{X')—fn{a)  = (x  — a)dn-, 
fn(a)  = F„(a)  -4-  g„(a)  = b- h gn{a)’, 


ïlff) 


Puis 

or 

donc 

(ii) 


O — a)  da=fn(x)  — b — gn(a), 

\x  — a\  dn £ |/,,(;r)  — 6|  + 1 ^„(a)|. 


Or  | gn  (a)  | ^ 8 d’après  ( 7 ) et,  puisque  l’on  considère 

comme  négligeable,  les  valeurs  de  x qui  vérifient  la  relation 

\f(x)~b\<$ 

vérifient  aussi  la  relation  | fn(x)  — b | < (3  et  alors  (11)  montre 
que 

| x — a | dn  < ^ -h  8. 

De  plus,  d’après  (2),  dn^.Dn — u \ donc 

(3-4-8 


(12) 


I x—  a | < 


Si  a est  une  valeur  approchée  de  a à a près, 

/ o \ 1 , | (3  -t-  8 

(i3)  \x  — a | < ha. 


— 4 — 

Remarques.  — i 0 Si  D„  — w,  on  s’en  tient  à n1  < x < rcw. 

2°  S’il  s’agit  de  trouver  une  limite  inférieure  et  une  limite 
supérieure  du  nombre  x défini  par  /(x)  = b , on  n’aura  qu’à  faire, 
dans  ce  qui  précède,  (3  = o. 


Si  maintenant  on  veut  tenir  compte  de  f (x) — fn(x),  posons 


04) 


/(*)  —fn(oo)  = R»(a?)i  * 


et  cherchons  une  limite  supérieure  de  | R„(:c)  |.  Pour  cela,  posons 
avec  M.  Peano  ( Formülario , p.  290)  : 

,®  (a,  b,  c)  = (b  — c)f(a)  -+-  (c,—  a)f  (b)  (a  — b)f(c)  (*), 

© ( C bj  c)  = (b  — c)f  ( t)  + (c  — t)f(b)  -1 -(t  — b)f  (c), 

4 O > = ( a — b)(a  — c)  y ( t,  b,  c)  — (t  — b)(t  — c)<p(a,6,c), 


a,  b,  c étant  des  nombres  donnés  quelconques  et  t une  variable. 
La  fonction  4 (t)  s’annule  identiquement  pour  t — a et  aussi  pour 
t — b et  pour  t = c,  car  cp  (6,  b,  c)  = o et  <p  (c,  b , c)  — o.  Donc, 
d’après  le  théorème  de  Rolle  (en  supposant,  pour  fixer  les  idées, 
a <<  b < c),  la  dérivée  <]/(£)  s’annule  pour  deux  valeurs  de  t com- 
prises l’une  entrera  et  b , l’autre  entre  b et  c;  donc  la  dérivée  , 
seconde  ( t)  s’annule  pour  une  valeur  tK  de  t comprise  entre  a 
et  c.  Or 


4"({)  ==  {a  — b)(a  — c)(b  — c)f(t)  — 2<p  ( a,  b,  c); 

donc 

(a  — b) (a  — c)(b  — c)  f" ) — îo(a,  b,  c)  = o. 


Le  raisonnement  et  la  formule  précédente  subsistent  quel  que 
soit  l’ordre  de  grandeur  de  a,  b:  c,  en  désignant  par  tK  un  nombre 
compris  entre  le  plus  grand  et  le 'plus  petit  des  trois  nombres  a, 
b , c.  Cette  formule  subsiste  encore  quand  deux  de  ces  trois 
nombres  sont  égaux,  car  alors  le  premier  membre  est  identi- 


(*)  C’est-à-dire 


?.(«,  à.  c)  — 


/(a)  a 
f(b)  b 
f(c)  c 
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quement  nul.  En  faisant  dans  cette  formule  a = x,  b n - J—  1 , 

c = ti  et  en  remarquant  que 

? (xi  n + i,  n)  = f(a ?)  — {x  — «)/(n  -f- 1)  — (/i  -h, i — a?)/( rc) 

= /0)  - /*/»(>)  d’après  (3)  = R*  (a?), 

on  a 

( 1 5 ) R« (a?)  = aï)  O - w — i)/r(*i), 

étant  compris  entre  le  plus  petit  et  le  plus  grand  des  trois 
nombres  .r,  n,  n- f-i.  Supposons#,  n et  n -f- 1 compris  dans  un 
certain  intervalle  (N,,  N2)  et  soit  m une  limite  supérieure  des 
valeurs  que  prend  j quand  L varie  de  N,  à N2  : 


(16) 

D’autre  part,  en  posant  X = n -f-  ^ -\-z, 

(*-»)(*!  «-0  = ^+1)  («-!)=*»- 1. 


3,  on  a 


, . . . i 3 

Etudions  la  variation  de  z2  — - quand  z croît  de  o à - : 


z 

17  3 

0 — 1 - 

210  2 

I 

1 , 3 

*2  — 

— - 0 0,2.4  - 2 

4 

4 4 

Or 

donc,  en  posant 
d7  ) 

(18)  si  \2\i~t 

(19)  si  M = « 


0,  4 < -,  ; 

4 


8W=f*’ 

1 < < 6 

ou  n — , 

5 D 


ou  /i , 


(20)  si 


z i - ou  « — 1 Sar.i/i  + 2. 


Gela  posé,  reprenons  les  deux  problèmes. 


| RnO)  | < fi, 

| R„0)|  < 3p, 
| Rre(a?)  | < 8 U. 
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Problème  direct.  — Si  n <.  x < n -f-  i et  qu’on  remplace  f(x) 
par  l’erreur  commise  est 

f(a ?)  — Fn(x)  = f(x)  — fn{x)  -»-/»( a?)  — F„(a?) 
ou 

O O’  f(x)—Fn{x)=Fn{x)  + gn(x). 

Donc,  d’après  (18)  et  (7), 


I/O)  — F„(ar)-|  <{*•-+- 8. 

Problème  inverse.  — b et  $ étant  des  nombres  donnés  (|S  > o), 
trouver  une  limite  inférieure  et  une  limite  supérieure  des 
nombres  x qui  vérifient  la  relation 

(22)  b — [3  </(a?)  < b -(-  (3  ou  | f{x)  — b | < 3. 

On  détermine,  comme  plus  haut,  un  nombre  n tel  que 

03)  V*^<V„+1 

■ n 

et  un  nombre  a tel  que 

(24)  F n(a)=b, 

d’où  n^a  <<  n- f-  1 . Mais  (21)  peut  s’écrire 


Or 

donc 


/O)  — b — gn(x)  — R„(a?)  = Fn(x)  — b . 


F„0)  — b = F„(a?)  — F„(a)  = (x  — a)  !)„; 


(25)  f{x)  — b — gn{x)'~  R n(x)  = (x  — a)  Dn, 

d’où 


(26)  \x  — a I ! fn(x)  — b | -h  | gn(x)  | -h  I R* O)  |. 

i°  Si 

Vn<(ù — [3  — S,  6 + ^ + ô < Vn+i, 

on  a 

V»-+-  S < b — p,  b -h  (3  < V„+1  — ; 


donc  n < x < n -)- 1 et  alors  de  (26)  on  tire,  d’après  (22),  (7) 


— 7 — 


et  ( 1 8 ), 
<*7) 

20  Si 

on  a 


\x  — a\  < 


0 -h  fx 


D, 


V,i  -1  < 6 — (3  — 8,  6-+-^-+-o<[Vn4-2, 

C/i— l r V,t-  1+  ô < b — (3,  è-hP<ÜV  /IH-2  Q v n- 1-2  i 


donc  n — j<^<n  + 2 et  alors  de  (26)  on  tire,  d’après  (22), 

(9)  et  (2°)> 

08)  1»-»|<P  + 3^8^ 

3°  Si,, en  appliquant  (28),  on  trouve  que  n — ^ <C  ^ 
on  pourra  tirer  de  (26),  d’après  (22),  (8)  et  (19), 


(29) 


. ^ S + 9.04-  3u 

D 


4°  Si  l’on  trouve,  en  appliquant  (28)  ou  (29),  que 


1 6 

n — - < a?  < /i  4-  - : 
5 5 


on  pourra  tirer  de  (26),  d’après  (22),  (8)  et  (18), 

2Ô  +'[JL 


(3o) 


a?  — a < 


J’ajoute  que  l’on  peut  appliquer  les  formules  (27),  (28),  (29), 
(3o)  sans  avoir  fait  de  constatations  préalables.  E11  effet,  la  for- 
mule (27),  par  exemple,  peut  s’écrire  x{  <^x  •<  a?2,  en  posant 

(3-f-Ô-hjJL  P-t-ÔH-W 

xx  — a — , x2  = a H-  i • 

Je  dis  que,  si  n <i  xK  <C  #2  <C  ^ -b  1 » cela  suffit  à prouver  que 
la  relation  |/(^)  — 6 | <;  (3  entraîne  x^<ix<^x2-  En  effet, 
d’après  (25), 

f{xi)  — b = (a?,  — a)  D„h-  &n(x\)  -+-  R«(j?i) 

= — P — & -4  R„(a:i). 


Or  n<^x{ 


par  hypothèse;  donc  gt)  (#,  ) b S d’après.  (7) 
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et  R/ï(#4  ) < p.  d’après  (18)  ; par  conséquent 

* b < — ^ — ô — p -t-  o H—  p 

ou /( x \ ) cb—  fi.  De  même,  b- bP</(#2)* 
Donc,  si  b — ^ < f{x)  < b -f-  [i,  on  a 

X\  < x < .r2. 


Calcul  de  /"O). 

Si  /(#)  — log  vulg.  on  a 

M 

avec  M = log  vulg.  e = o,43 . 


donc 


\f'(x)  | < 5. 10  7 pour  x > 1000. 


Pour  les  tables  des  log  tang  des  arcs  de  centigrade  en  centi- 
grade, 

Y HZ  OC 

f(x)  = log  vulg.  tang^-  avec  y=— -, 

Mit2  cosv 


,,  = » 

J X ^ 104  siny 


f"(x)  - 


108  sin2^ 

Si  600  < X <9400,  [/"(a:)}  <12.10-7; 
Si  3ooo  < x < 7000,  \f"(x)  | < 4.  io~8. 


On  peut  aussi  calculer  /"(.r)  au  moyen  des  tables  elles-mêmes 
quand  f"{x)  varie  constamment  dans  le  même  sens  dans  l’in- 
tervalle  considéré.  En  effet,  posons 

d(x)  = f(X  H-  1)  —f(x). 

On  a,  d’après  la  formule  des  accroissements  finis, 

d{n-+-p)  — d(n)  = p d\xi)  = />[/' Oi -h  1)  — f {x ,)]  = pf'(x2) 
avec 

n < X\  < n x 1 < x 2 < x 1 -t-  1 . 

Donc 

dn+p  — dn=  pf"(x 2),  n < x2  < n —1—  p 1 . 

De  même, 

du- , — = pf{xA),  n — p — 1 < x3<  n. 
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Donc,  si  f"{x)  va  constamment  en  croissant, 

dn— i — dn-p—\  < pf"(n)  <C  dn+p  dn 

ou,  d’après  (2), 

D„_,  — Dn-p-,  — iu  <pf"(  n)  < Drt+jo—  D„+2  u. 

Par  exemple,  pour  les  tables  à cinq  décimales  des  log  tang  des 
arcs  de  centigrade  en  centigrade, 

D499  = 87  a,  D599  = 73  U r=  D6oo,  D700  = 63  U 

avec  u = io~5.  Donc 

— 16. 10  7 </"(6oo)  < — 8.  i o-7. 

Le  lecteur. — Quand  j’ai  passé  le  baccalauréat,  on  m’a  demandé 
de  faire  un  calcul  de  logarithmes,  avec  toute  la  précision  que 
comportent  les  tables  : que  voulez-vous  de  plus  ? 

L’auteur.  — Pour  passer  le  baccalauréat,  il  suffit  d’observer 
les  règles  du  jeu  du  baccalauréat;  mais,  quand  tu  voudras  faire 
un  calcul  sérieux  avec  des  tables,  tu  seras  bien  obligé  de  te  rendre 
compte  du  degré  de  précision  que  comportent  ces  tables. 

Le  lecteuh.  — Vos  formules  sont  bien  compliquées;  je 
préfère  la  formule  que  donnent  les  bons  auteurs. 
r 

L’auteur.  — Il  est  difficile  de  trouver  des  règles  qui  soient  à 
la  fois  simples,  exactes  et  générales.  On  peut  essayer  de  faire  des 
règles  particulières.  Ainsi,  pour  f(: r)  = logvulg.#, 

M 

/O)  — fia)  = 0 — a)  — , 

X\ 

x x étant  compris  entre  a et  x. 

En  posant,  comme  plus  haut,  F n(a)  = b,  on  en  tire 

x — a f(x ) — b -h  F n(a)—fn(a)  -*-/„( a)  — f(a) 

xl  M 

Or  = 2,3°  . . . < 2,5  et  |Fw(a) —/„(«)!<  5 d’après  (7); 

donc,  si  \f(x)  — b | < et  que  /»(«)  — /(a)  soit  négligeable, 
on  a,  en  désignant  par  x ' une  limite  supérieure  de  a et  de  a?  et, 


par  suite,  de  , 


En  mettant  x au  lieu  de  x\  3 au  lieu  de  2,5,  et  en  négligeant  S, 
tu  pourras  dire  que  V erreur  relative  commise  sur  un  nombre 
calculé  au  moyen  de  son  logarithme  est  moindre  que  le  triple 
de  V erreur  absolue  commise  sur  ce  logarithme. 

Le  lecteur.  — Pourquoi  n’avez-vous  pas  parlé  de  la  différence 
seconde,  ni  de  la  formule  classique  de  JNewton? 

L’auteur.  — En  désignant  la  différence  seconde  du+t  — dn 
para^2),  la  formule  de  Newton  s’écrit 

f(x)  = f{n)  -h  (x  — n)dn- f-  - (x  — n)  (x — n — i 

v 2 

Pour  se  servir  de  cette  formule,  il  faut  la  limiter  à un  certain 
terme  et  ajouter  un  terme  correctif.  En  s’arrêtant  au  deuxième 
terme  et  en  désignant  le  terme  correctif  par  Rn(#),  on  a juste- 
ment la  formule  (i4)* 
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